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Abstract. We investigate the viability of having dark matter in the minimal left-right 
symmetric theory. We find the lightest right-handed neutrino with a mass around keV as 
the only viable candidate consistent with a TeV scale of left-right symmetry. In order to 
account for the correct relic density with such low scales, the thermal overproduction of the 
dark matter in the early universe is compensated by a sufficient late entropy production 
due to late decay of heavier right-handed neutrinos. We point out that the presence of the 
right-handed charge-current interactions, operative around the QCD phase transition, has a 
crucial impact on the amount of dilution, as does the nature of the phase transition itself. 
A careful numerical study, employing the Boltzmann equations, reveals the existence of a 
narrow window for the right-handed gauge boson mass, possibly within the reach of LHC 
(in disagreement with a previous study). We also elaborate on a variety of astrophysical, 
cosmological and low energy constraints on this scenario. 
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1 Introduction 

The Left-Right (LR) symmetric theories [1-4] have over the years emerged as one of the main 
candidates for the theory beyond the Standard Model (SM). Its great achievement is a pre- 
diction of non- vanishing neutrino mass, whose smallness gets naturally tied to the maximality 
of parity violation of weak interactions, through the seesaw mechanism [5-9]. This model 
offers spectacular signatures at colliders such as the Large Hadron Collider (LHC): direct 
lepton number violation through the production of the heavy right-handed (RH) neutrinos 
and its subsequent decay, giving a final state with same-sign dileptons [10], the high-energy 
counterpart of the neutrinoless double beta decay [11, 12] and lepton number violating de- 
cays of light mesons [13]. Moreover, the Majorana nature of the RH neutrino manifests itself 
through the equal production of charged leptons and its anti-particles. A confirmation of the 
claimed observation of the neutrinoless double beta decay, could even require the LR scale 
to be tantalizingly close to the LHC reach [14, 15] (for a review, see [16, 17]), if cosmology 
constraints keep pushing down the sum of neutrino masses [18-20] and thus disfavoring the 
contribution due to neutrino mass [21]. 
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There is circumstantial evidence that Dark Matter (DM), which composes about eighty 
percent of matter energy density in the universe, is in the form of particles. As is well known, 
the minimal SM fails to provide a DM candidate. At first glance, LR theories could do the job, 
for they introduce new, potentially stable particles. In particular, the lightest RH neutrino, 
if light enough, could easily be cosmologically stable, since its Yukawa couplings to the light 
neutrinos can be arbitrarily small. In a sense, the LR theory can be viewed as a natural 
framework of DM. Notice there is no need for artificially imposing any Z2 symmetry at high 
energy - the SU{2)ii gauge interactions would simply break it. Such gauge interactions 
cease to be operative in the decay of the lightest RH neutrino, if it is the lightest fermion 
that couples to Wr. An approximate Z2 symmetry emerges at low energy, if all its Yukawa 
couplings are negligibly small. 

A warm DM candidate [22-25] with a mass around keV, while it works as well as cold 
DM for the large scale structure formation, can suppress the structures on smaller scales via 
free streaming [26]. This scenario is of particular interest as a solution to the problems of 
very cuspy halo profiles and over-populated low-mass satellite galaxies, usually predicted by 
cold DM. The idea of having RH neutrino as warm DM candidate with a mass around a keV 
was introduced around thirty years ago [25, 27]. Due to the presence of gauge interactions, 
one expects the RH neutrino playing the role of DM to have a similar relic number density 
as the one of the light neutrinos, if the scale of LR symmetry is not far above the electroweak 
scale. In [25], the first cosmological bounds on stable heavy neutrinos, charged under a 
new gauge symmetry, were studied and the problem of their potential over- abundance was 
stressed. Ref. [28] offered a nice way out, by today a text-book scenario [29], through the late 
entropy injection due to the decay of a heavier long-lived particle, for example the heavier 
RH neutrino [30]. 

Using the idea of [30], a few years ago Bezrukov et al [31] performed studied this issue 
in the context of LR theories and argued that one cannot obtain the correct DM abundance 
unless the mass of the RH charged gauge boson Wr is above 10-16 TeV, far from the LHC 
reach. This unfortunate result made it irresistible for us to reconsider their analysis with 
great care. While we agree with the basic mechanism presented in [31], our analysis reveals 
an additional window for the Wr mass around roughly 5 TeV, possibly within the LHC reach. 
This is the main result of our paper, whose importance cannot be over-emphasized. 

The key point in our work of realizing such low scale LR symmetry is to take advantage 
of the QCD phase transition, where the number of relativistic degrees of freedom changes 
dramatically. Depending on the flavor structure of the their gauge couplings, the RH neu- 
trinos decouple at different temperatures, which could be separated enough to lie before and 
after the transition. We show that this plays a crucial role in producing large enough amount 
of entropy, in order to dilute the DM relic abundance towards the acceptable range. 

In the decoupling limit, when the LR scale is large, our picture crosses over smoothly to 
the so-called i/MSM scenario, studied extensively over the years [30, 32, 33] (for a review on 
the topic of light sterile neutrinos, see [34] and references therein). In the z^MSM case, one has 
only the SM augmented with RH neutrinos and Yukawa interactions take over the role of the 
RH gauge interactions, leading to the sterile neutrino picture. This transition is quantified 
carefully in Sec. 3.6. In other words, the LR theory cannot fail to account for dark matter, 
if its scale is high enough. However, the theory then stops being directly verifiable and loses 
most of its phenomenological appeal. This is why we find the existence of the narrow band 
for not-so-heavy Wr important enough to warrant another paper on the subject. 

In order to ease our reader's pain, we have decided to make this presentation as ped- 
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agogical as possible and therefore we wih start from scratch in presenting our work. In the 
fohowing Sec. 2, we will review the essential features of the LR symmetric model and sum- 
marize the present-day theoretical and experimental status. Section 3 presents all the central 
ideas behind this works in an accessible way, using rough estimates, which paves the way for 
Sec. 4, where a detailed numerical analysis is carried out. In Sec. 5 we discuss additional 
astrophysical and cosmological constraints, together with the limit coming from the search 
for neutrinoless double beta decay. We conclude in Sec. 6. 



2 Minimal Left-Right Model: A Telegraphic Review 

The left-right symmetric model is based on the SU(2)l x SU{2)b, x U{1)b~l gauge group 
(suppressing color), supplemented by a symmetry between the left and the right sector [1-4]. 
Quarks and leptons come in symmetric representations 

The Higgs sector of the minimal model [5, 6] consists of a bidoublet $ = (2/,, 2^2, O^-l) and 
two triplets, = (3l, Ir, 2b^l) and Ar = (1/,, 3r, 2b-l) 

Hereafter, we refer to this setup as the minimal LR standard model (LRSM). The symmetry 
breaking in the model [35] is characterized by the following vacuum expectation values (vev) 
($) = diag {vi, V2), (A^ ^) = vl^r, which have a hierarchical order vf^ <^ v"^ = vf + V2 <^ vj^, 
with V = 245 GeV. For a recent summary of the main features, see [36-38]. The resulting 
masses of the heavy gauge bosons are Mw^ = QVr and 

Mz^^-V^Mw^. (2.3) 

Notice that Zlr is appreciably heavier than Wr^ which turns out to be important in the 
study of the RH neutrino freeze-out. The vev's of the bi-doublet give masses to the charged 
fermions and Dirac masses to neutrinos. The triplet vev vr gives directly Majorana mass to 
RH neutrinos, which results in the type-I seesaw mechanism [5-9], while vl independently 
gives Majorana light neutrino masses, as in the so-called type-II seesaw [35, 39, 40]. As we 
show in what follows, the seesaw picture is a must for the DM scenario to work. 

The charged and neutral currents relevant for our study, up to tiny v^/vl^ < 10^'^ 
corrections (see below for the limits on the LR scale), are 



^'r V^/ r. 



+ h.c. , (2.4) 



Cnc = -^=L==Z^^J^, [T.R + tan2 9w{nL - Q)] f + :^A%NN (2.5) 
V 1 - tan^ 9w ^^^iwji 

where A^'s are defined as RH neutrinos mass eigenstate states, and we have taken Qr = Ql = 9, 
appropriately for a LR symmetric theory. The mixing matrices and are the right- 
handed analogues of the left-handed PMNS and the CKM mixing matrices with elements 
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Particle 


Final state 


Lower limit 


Collaboration 


Comments 


Wr 
Wr 
Wr 

Zlr 
Zlr 

Al 
Ai 


jj 

e/ii + N 
Ujj 


L5 TeV 
2.5 TeV 
< 2.5 TeV 
~ 2 TeV 
~ 3 TeV 
100-355 GeV 
70-90 GeV 
45 GeV 
113-251 GeV 


CMS [44] 
CMS [45] 
ATLAS, CMS [46, 47] 
ATLAS [49] 

LEP [51] 
ATLAS [54] 
LEP [58] 
LEP [51] 
ATLAS [54], CDF [56] 


independent on N mass 
light (missing energy) 
heavy Major ana N [48] 
see [50] 
indirect, see [52, 53] 
spectrum dependent [55] 
chargino search [57] 
Z-boson width 
flavor dependent 



Table 1. A summary of limits on the mass scales of the particles in LRSM from collider searches. 



Vlfi and V^,. Of course, there is no reason that gR = ql relation should hold exactly at the 
scale of interest; there could be easily a variation if the LR symmetry is broken at a high 
scale. The small effects due to renormalization group equation running would change none 
of our conclusions. 

Notice that Zlr, besides being heavier than Wr, has also smaller couplings to A^'s. 
Thus, a RH neutrino that couples only to Zlr will decouple earlier in the thermal history 
of the universe for a given LR scale. Since a warm DM candidate in this kind of a setup is 
typically overproduced, it will turn out desirable to profit from this fact and decouple it from 
the Wr at relevant temperatures. 

The most stringent theoretical limit on the LR scale is derived from neutral kaon mix- 
ing [41, 42], and the latest studies set a lower bound M\Yji > 2.5-4 TeV [36, 37, 43], de- 
pending on the choice of the LR symmetry, charge conjugation or parity, respectively. The 
experiment, however, is now catching up and the theoretical constraints are becoming obso- 
lete. Direct searches are continuously pushing up the limits on mass scales in the LRSM and 
we summarize them in Table 1. The window around ~ 5 TeV that will emerge from our DM 
study is comfortably above all the current theoretical and experimental bounds. It is also 
noteworthy that the second Higgs doublet belonging to orthogonal to the SM-like one, 
must be heavier than ~ 10 TeV due to the contribution to tree-level flavor changing processes. 
For a recent complete study of a variety of flavor processes in the LRSM, see [38]. 

Moreover, Table 1 tells us that most of the states have masses around or above the 
weak scale. The only exception are the RH neutrinos N and the neutral Higgs A^j, which 
behave like singlets under the SM gauge group and are not very much constrained by collider 
searches. They may be long-lived and are therefore potential DM candidates. In the next 
section we study which, if any, can actually do the job. 

3 A Tale of Three Right-handed Neutrinos 

We discuss here the history and role of RH neutrino as DM in the early universe, in the 
context of TeV scale LRSM. It contains the essence of what is going on: the lightest is 
presumably the DM due to its longevity, while the heavier one(s) should make sure that its 
abundance is correctly accounted for. 

This section is the core of our work; it is here that our reader will find the central ideas, 
albeit simplified. Thus, we urge her to postpone the coffee break until having gone through 
it. The technicalities required for a precise quantitative picture are left for the Sec. 4. 
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3.1 Warm Dark Matter Candidate 

Let us now go through the hst of potential candidates for the dark matter in the LRSM. From 
the discussion in the previous section, it is clear that the neutral components belonging to 
the Higgs bi-doublet $ and the triplet A/, are heavy enough to decay on collider time scales, 
and are thus ruled out from the start. 

• One possibility is the neutral component of the SU{2)pi triplet, ReA^j. Being a SM 
singlet, it is allowed to be as light as one wishes, and if it is lighter then the RH 
neutrinos, it will decay into two photons, with a rate 

r...„ - ? (^V - 10-» GeV C^)' ('^^)\ (3.1) 

^R^^^ Svr V47r/ \Mwn) \/2 VkeV/ V ^%if, J 

where 49 is the loop function squared. This approximate formula takes into account 
only the dominant contribution due to the heavy charged Wr. This vector boson 
dominance resembles the situation in the SM, where the gauge contribution is roughly 
an order of magnitude bigger than the rest, i.e. the fermionic one. Here, there is no 
fermionic contribution, but instead the one from the charged scalars, which is typically 
much less than the gauge boson one. 

There is a stringent lower limit on the stability of a radiatively decaying dark matter 
particle [59], ta > 10^^ sec, or equivalently T < 10"^'' GeV, up to uncertainties in 
astrophysical parameters. This pushes the scale of LR symmetry far away from the 
LHC hope into the despair of no direct detection. It is somewhat surprising though, 
that the A^j could be a viable dark matter candidate in an 50(10) grand unified theory, 
where the LR symmetry breaking scale lies preferably around 10^^ — 10^^ GeV [60, 
61] (admittedly it would have to be incredibly light, creating yet another hierarchy 
problem). In any case, this possibility would take us far away from our search for a low 
scale LR symmetry. 

• Thus, for low LR scale, the only viable candidate left is the lightest RH neutrino, to 
which for definiteness we refer as A'^i. For its mass below the pion mass, the decay 
channel mediated by Wr closes quickly. It can only be destabilized by Dirac Yukawa 
couplings due to the mixing with left-handed neutrinos. This mixing leads to its decay 
to a light left-handed neutrino and a monochromatic photon, or three light neutrinos. 
Again, the X-ray constraints implies such mixing to be tiny [62, 63] 

0?<(1.8-3.1)xlO-5f— v. (3.2) 

In the rest of this work, we thus study the exciting possibility of A^i playing the role 
of dark matter. The most reliable cosmological lower limit on the DM mass is derived by 
considering the phase space density of compact objects, which is around ~keV scale [64- 
66]. For M\Y^ lying in 1 — 10 TeV region, one may worry whether such small mass of A^i is 
consistent with the seesaw formula, with radiative corrections included. The real question is 
the radiative stability of the neutrino Dirac mass, studied in Ref. [67], with a conclusion that 
it can be as small as a few eV. This in turn implies (via the seesaw) that the RH neutrino 
is allowed to be as light as 10 — 100 eV. In other words, a keV RH neutrino is a perfectly 
natural choice for a DM candidate. 
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The main obstacle we have to face when M^i/^ lies in the TeV region, turns out to be 
the over-abundance of A'^i [25, 31], because the SU{2)ji gauge interactions keep it in thermal 
equilibrium when the temperature is high. Intuitively, one expects A'"i to decouple at a 
temperature not far above the usual decoupling temperature of light neutrinos, so that its 
density is also similar to that of light neutrinos - a disaster for a particle with a mass above 
keV. Obviously, the heavier it is, the more severe a problem this becomes. The bottom line 
is that one is pushed to the picture [27] of warm dark matter, with a mass in the keV range. 
This is welcome due to the need to suppress the small-scale structures. 

Furthermore, it has been pointed out that supernovae cooling imposes a tight constraint 
on the flavor structure of weakly interacting species lighter than lOMeV [68]. In the context 



of LRSM, this implies M\Yj^ > y | Vj}\ x 23 TeV [69] , which means that the electron component 
of A^i in the RH charged current has to be roughly below 1%, if LR symmetry is close to the 
TeV scale. 

In what follows, we shall stick to the keV warm dark matter RH neutrino and pursue 
its implications to the bittersweet end. 

3.2 Thermal production via freeze out 

The presence of new gauge interactions in the LRSM has a major impact on the thermal pro- 
duction of dark matter in the form of RH neutrinos. If the universe starts from a sufficiently 
high temperature, these interactions will keep them in thermal equilibrium via scatterings 
with the SM fermions. On the other hand, it turns out that the Dirac Yukawa couplings 
of a keV RH neutrino are never large enough to bring it into equilibrium above the elec- 
troweak scale. Also, after the electroweak symmetry breaking, the mixing between A'^i and 
light neutrinos in matter are suppressed by the finite density potential [70, 71], compared to 
the vacuum mixing angle 9i. Therefore, the SM weak interactions cannot bring RH neutri- 
nos in equilibrium [63]. Although the direct thermal production via the mixing is negligible, 
non-thermal contributions due to accumulative oscillations [27, 72-74] may become signifi- 
cant, as in the case of z^MSM. However, this effect is sub-dominant compared to the thermal 
(over-) production in the LRSM as discussed below. 

Let us quantify more precisely how the dark matter RH neutrino A'^i with TeV scale 
gauge interactions typically gets over-produced in the usual thermal history. The freeze-out 
temperature Tf can be estimated by the out-of-equilibrium condition T = H, where T is the 
annihilation rate of A^i and H is the Hubble parameter. The interactions that keep them in 
thermal equilibrium are scatterings of with the light SM fermions, mediated by the heavy 
gauge bosons, Wr and Z^r in the LRSM. In the radiation dominated era, we have [75] 



Generally, for any RH neutrino that decouples from equilibrium while still relativistic, one 
can obtain the freeze-out temperature Tf as a function of Mpv^ and the number of degrees 
of freedom g^^iTf) at that time. 

For the Wr in the TeV region, the freeze-out temperature is around 





(3.3) 




(3.4) 
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where the normahzation reflects the fact (7=,, (400 MeV) ~ 70. For a RH neutrino A^, which 
freezes out while relativistic, the number per entropy density is approximately 



UN 135 C(3) 



(3.5) 



The yield Yjsf, being thermally conserved quantity, turns out to be useful in what follows. 
It may be suppressed if Tf is higher, which can happen if neutral current interactions play 
the dominant role in the freeze-out, as pointed out above. Such a situation occurs, when 
the presence of the charged lepton corresponding to a RH neutrino in the thermal bath is 
Boltzmann suppressed. With the above values of the freeze-out temperature, the optimal 
option is to couple the DM candidate A^i predominantly to r. This is the first crucial 
ingredient towards determining the flavor structure in the RH charged currents necessary to 
accommodate DM in this theory. Remarkably enough, when the dust settles we will end up 
with a completely determined flavor structure of the charged current. 

In any case, a keV A'^i is still relativistic at such temperatures, therefore its relic abun- 
dance today is roughly 



where we have used today's entropy density s = 2889.2 cm~^, critical density pc = 1.05368 x 
10~^/i^ GeV/cm^ and h = 0.7. This is to be contrasted with the observed dark matter relic 
abundance of the universe [83] 



at 3cr confidence level. Clearly, the estimated contribution of A'^i in Eq. (3.6) over-closes the 
universe, when Wr lies in the TeV region, by at least a factor of ~ 12.5 x (m^r^/l keV). Since 
there is no room in the minimal LRSM for ~1000 degrees of freedom in g^[Tfi), this problem 
cannot be solved by simply raising the scale Mpv^j. Such a possibility of a huge number of 
ad-hoc new states is rather unappealing in any case. 

3.3 Late entropy production 

The only way out of this impasse is to dilute the number density of A'^i by entropy production 
due to the late decay of some massive particle which dominates the universe [28] . Such a late 
decay should inject relativistic light SM particles that quickly equilibrate with the thermal 
plasma and "reheat" the photon temperature. In turn, it takes longer for the photons to cool 
down to present-day temperature and the number density of DM is effectively reduced. In 
order for the dilution to work, the temperature of A'^i should not increase, and it is therefore 
crucial that A''! itself is not a decay product of the heavy decaying particle. 

In order to release a substantial amount of entropy, such a particle is required to be 
long-lived, with lifetimes up to a second. In the minimal LRSM, the only particles relevant 
for the late decay are the following. 

• As discussed at the beginning of Sec. 3.1, the neutral component of the SU{2)ji triplet 
ReA^j can be as light as one wishes. For TeV scale LR symmetry, it has to be lighter 
than MeV in order to live as long as one second. This makes it too light to play any 
significant role in entropy production. 




(3.6) 



n^M = 0.228 ± 0.039 



(3.7) 
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• The heavier RH neutrinos are the only remaining viable candidates and they play 
the role of diluters in this scenario [30]. This section focuses on determining their 
characteristics, required for adequate relic abundance of DM. Here, we refer to the 
diluter as a generic N, and use explicit indices (2,3) when necessary to specify the 
flavor structure. The DM RH neutrino is always A'^i. 

In order to achieve a sufficient dilution, the mass of the diluter m^v should not exceed its 
freeze-out temperature Tj, otherwise the yield in Eq. (3.5) receives an additional Boltzmann 
suppression factor g""*^/-^/. If the gauge interactions of RH neutrinos are universal, their 
freeze-out temperatures is of the same order (we return to this point in Sec. 3.5) and similar 
yields for all RH neutrinos Yat^ are expected. 

As the temperature of the universe drops, a sufficiently massive and long-lived RH 
neutrino can temporarily dominate the total energy density. After decays, the energy 
density is transferred into that of radiation. In the sudden decay approximation, all the A^'s 
decay at t ~ rjv and "reheat" the universe to the temperature Tr [28] 

/I \ 1/2 

Tr ~ 0.78g^{Tr)-^/^y/TNMp ~ 1.22MeV ( ^ j • (3-8) 

In order to start the Big Bang Nucleosynthesis (BBN) with a correct proton-neutron number 
ratio, Tr should be larger than about MeV, which gives an upper bound on diluters' lifetime 
Tjv < 0(1) second. 

Using energy conservation m]\rniy{T]y) = miyYiys = PR{Tr), or in the other words 

3 

niNYjqSbefore = -^SafterTr , (3.9) 

the dilution factor, defined as the ratio of entropy before and after the decay, in the sudden 
decay approximation (no volume change) becomes roughly 

S ^ ^ ^l.8{g,{Tr)fl^ , (3.10) 

'-'before ■Sbefore ^jY ]\jM^ 

where Fat = r^^ and Mp = 1.2 x 10^^ GeV is the Planck scale. If such entropy production 
happens well after A"! froze-out, the relic density calculated in Eq. (3.6) will be reduced by 
the dilution factor, VLn^ — )• O^Vi = ^N]_/S. For a reheating temperature Tr around MeV, 

A . , / mAT, \ /1.85GeV\ /lsec\^/V5*(r/-2 3)\ 

«„. .(0,228 + 0^039) (^)(^^)( — ) [^^)- (3,11) 

The dilution factor in Eq. (3.10) is proportional to the lifetime of the diluter A^, which 
should be long enough to make a sufficient impact on the relic density of A^i. On the other 
hand, a late decaying particle with a lifetime longer than about a second may threaten the 
success of BBN. Therefore, the optimized situation is to have r^r ~ 1 sec (or Tr ~ MeV, as 
assumed above) which, depending on the scale of LR symmetry, narrows down the mass 
range of the diluter A^, since the same gauge interactions that govern the freeze-out of RH 
neutrinos are responsible for the decay of A^. 

Let's recapitulate once again the logic of our search for the light Wr accessible to 
experiment. This immediately narrows down the freeze-out temperature of A^'s in the few 
100 MeV range, as seen in Eq. (3.4). In order to avoid the Boltzmann suppression, the 
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diluters mass should be below Tf, which in turn brings back the problem of over- abundance 
of DM, since the third term in Eq. (3.11) becomes large. We are back to square one, it seems. 
However, there is a potential way out; make the last term small by separating the freeze-out 
temperatures of the diluters and the dilutee. It turns out that the nature of the QCD phase 
transition plays an essential role in this, as discussed in the coming section. 

A lifetime of around one second restricts the mass of to a narrow region in the few 
100 MeV range for a Wr in the few to ten TeV region. Depending on its mass, N decays 
via heavy Wr either predominantly mainly into a lepton plus two light quarks (which later 
hadronize), with a lifetime 

r(N- m - — ^ - Isec f f .3 12) 

or (if rriN ^ ttIt^ + rrii, i.e., near the threshold) into a lepton and a pion, with a lifetime [15] 



1 sec 



mN ( Mw^ y ( 0.002 

250MeV/ V^TeVy \f{x^,x^ 



1 /2 

where f{x^, x^) = [(1 - xjf - xl{l + x})] [(l - {x^, + x^f) (l - (x^ - Xif)] and x^^^. = 
TriT^^i/mN ^. In the above estimates, we took V^^ ~ ^ud^^ ~ ^- Fi'o™ the above expressions 
for N lifetime, it is clear that the best bet to have M^^^ in the TeV region is to have the 
pionic decay dominant, together with a final-state phase space suppression. This narrows 
down niN to lie around 

rriN ~ "T-TT + rni , (3.14) 
together with a lepton mixing V^^ ~ V^^ f« 1, up to ~ 1% (see Fig. 3). 

3.4 Fixing the flavor structure 

Let us take a closer look at the leptonic flavor of the diluting RH neutrino N in Eq. (3.14). 
In order to successfully dilute the DM relic density in Eq. (3.11), one would naively conclude 
that £ = r is favored, since the mass of the diluter is in the right ballpark. However, there 
are several serious drawbacks related to this channel. 

First and foremost, in order to avoid the non-relativistic Boltzmann suppression, the 
freeze-out temperature of N should be bigger than its mass. In order to lift this suppression 
a high freeze-out temperature is needed, which which for m^r around 1.85 GeV, requires the 
Wr boson mass to be heavier than about 15 TeV [31], see Eq. (3.4) (see also Fig. 3). Being 
far out of the LHC reach, this is outside of the region of our interest. This by itself prevents 
any diluters to couple predominantly to the tau. 

There is yet another reason that strengthens this result. Namely, the decaying 
might have an appreciable branching ratio to A^i, which makes the task of dilution more 
challenging. Namely, the cosmological lower limit on the warm DM mass depends on its 
free-streaming length, which is proportional to the average energy {pni) at injection [84]: 
\fs — IMpc ((pAr-^)/(p,y)) (keV/m^Vi) 5"-*^/^. The analysis of the Lyman-a forest demands 
the DM free-streaming length not much longer than a Mpc. Unless the mass of A^ happens 



^Due to a small mixing angle 9% ~ m^/raN, the A'' — > channel turns out to be subdominant. 
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to lie within few MeV above the production threshold, the A^i in the decay product is always 
too energetic. The free streaming of A^i will erase any structure at scales larger than a Mpc, 
unless misf-^ is much bigger than keV. From Eq. (3.11) such a heavy mass is disastrous for 
relic density, for a TeV LR scale. 

Therefore, in what follows, we focus either on the electron or the muon flavor channel for 
the diluters. In this case, the decays of A''2,3 to (an energetic) A^i are kinematically forbidden 
if A'^i couples predominantly to r, which is the preferred freeze-out scenario anyway. Namely, 
in that case the freeze-out temperature of A'^i is dictated by the new neutral currents only, 
suppressed compared to the charged ones, which in turn makes A^i decouple before the 
diluters. At this point, a clear picture of both the flavor structure and the mass spectra of 
RH neutrinos emerges 

/O 1\ niN, keV, 

Vf« 010 , mN,-m^ + m^, (3.15) 

\1 0/ TTlATg m,r + me. 

Notice that any departure from diagonallity in the 2 — 3 sector effectively reduces the masses 
of diluters and therefore decreases the dilution factor. In other words, we end up with a 
flavor-diagonal RH leptonic mixing matrix, to which we stick for the rest of this paper. 

Such a flavor structure is completely different from the PMNS mixing matrix in the left- 
handed sector. In the case of the type II seesaw, due to the LR symmetry, the flavor structure 
in the left and right sectors has to be identical. This means clearly that the contribution 
from the type II seesaw can only play a subdominant role. Let us see now how the type-I 
seesaw dominance for neutrino masses ends up being consistent with the diluters' lifetimes. 
In fact, the partial lifetime of a RH neutrino induced by the Dirac yukawa coupling in the 
type I case is 

r,= (^y(^yr,^(^y^r,^( ^li^y 10^ sec. (3.16) 



For the diluters with mass around 0.1 GeV, this is certainly a subdominant channel in the 
presence of TeV gauge interactions and they can fully participate in the seesaw. This is in 
contrast to the case of mjv > GeV, where the Dirac mass has to be suppressed to ensure 
enough dilution, in which case one would end up with type II dominance (as in Ref. [31]). 
Whereas in the latter case, the W^^ is necessarily heavy, far above the LHC reach, the 
appealing aspect of the former case is an emergence of a light window, on which we 
focus. 

The bottom line then is, we end up with fairly light RH neutrinos and negligible leptonic 
mixing, which eliminates completely both the lepton number violation at colliders and lepton 
flavor violation. A positive future result in any such process would kill the picture envisioned 
here, and dark matter would have to come from elsewhere. 

3.5 A Window for low scale LRSM 

For obvious reasons, we here focus on the attractive possibility of having VFr as light as 
possible. From Eq. (3.11), it is clear that the ratio of g^{Tf2^3) / g*{Tfi) should be small 
enough in order to compensate for the smallness of niN. This amounts to decoupling A^i 
as early as possible, compared to N. In the LRSM, the main processes that keep any RH 
neutrino in thermal equilibrium are 

Nr^ud, Nu^i+d, Nd^£+u, (3.17) 
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Figure 1. Left. Thermally averaged reaction rates for the processes that dominate the decoupling 
of A^i (red) and N2 (blue), for parameters = 5 TeV, = 0.25 GeV. Also shown is the Hubble 

expansion rate multiplied by the thermal number density of N2 (black curve). Right. A sharp change 
in the evolution of around the QCD phase transition temperature. The dot-dashed (dashed) curve 
corresponds Tqcd = 400(150) MeV with a second-order QCD phase transition (taken from [85]), while 
the solid line is an interpolation in between with Tqcd = 350 MeV and the transition close to first 
order. The red and blue points are the freeze-out temperatures of iVi and N2, respectively. 

mediated by Wn, and 



mediated by both Wr and Zlr. The relevant cross-sections and their thermally averaged 
reaction rates are given in the Appendix B. We find that the single- annihilation rate 
dominates the pair annihilation at least by one order of magnitude. The main reason is the 
suppression of the quark couplings to Zlr, as mentioned above. 

A crucial point to note is that, in the single-A^ annihilation processes, when the tem- 
perature drops below the mass of the corresponding charged lepton, the reaction rate starts 
to be Boltzmann suppressed, by more than one order of magnitude for temperatures around 
and below the QCD phase transition. ^ Decoupling A'^i earher than N2 therefore requires a 
particular flavor structure of the right-handed leptonic mixing matrix V^. Namely, the A^i 
should couple predominantly to the heaviest charged lepton, r, i.e. ~ 1 (cf. Eq. (3.15)). 
In such a case, when the temperature drops below the r mass, A'^i is maintained in the equi- 
librium only through the neutral current interactions, governed by Zlr. On the other hand, 
A^2,3 now couples to or e and keeps annihilating through the charged current interactions, 
too, and thus decouples later than A^i. In practice, the difference between the two decoupling 
temperatures can be as large as a few hundred MeV. 

Such a difference between the temperatures turns out to be important if the thermo- 
dynamical nature of the universe changes at the same time. For Wr around a few TeV, the 
freeze-out temperatures of A'^ lie around a few 100 MeV. This is precisely in the vicinity of the 
QCD phase transition, where a sharp drop in occurs, as most hadron states are becoming 
non-relativistic (right panel of Fig. 1). In case the total entropy is conserved during the 
transition, the temperature of the thermal plasma gets "reheated" due to the change of g^f. 

^There is a subtlety in the treatment of relevant annihilations at temperature close to the QCD phase 
transition. The issue is whether one should work with quarks or light mesons. Admittedly, the quark picture 
we opted for brings in uncertainties that cannot be easily quantified. In practice, since we ask the RH neutrinos 
to freeze out before or immediately after the transition, we believe such uncertainty is small. 




(3.18) 
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Figure 2. The yield of N\{N-2) at freeze-out is shown in magenta (blue), depending on the mass of 
Mv/ji- Different types of lines correspond to different QCD temperatures as denoted in the leggend. 
Notice that the nature of the QCD phase transition becomes irrelevant once Mvi/^ > lOTeV, since 
the freeze-out temperature is above GeV. See also the right panel of Fig. 4 for more precise such 
dependence evaluated using Boltzmann equations . 



For the sake of illustration, we show first what happens in a simple example with a 
single diluter, which we choose to be A''2 due to its larger mass and a bigger impact on the 
dilution factor. If the left-right scale is chosen, such that A'^i freezes-out before and A''2 after 
the phase transition (see Fig. 1), will feel this reheating which will enhance its yield. The 
relative number density between two relativistic RH neutrinos is 

^ ~ ^-(^^^^ (3 19) 

and can be as large as 3 — 4 for M^y^ < 10 TeV, as shown in Fig. 2. This will eventually result 
in more entropy release during the later decay of A''2, and it is manifest in the reduction of 
the ratio g*{Tf2) / g*{Tfi) in Eq. (3.11). Although in this simplified picture the enhancement 
factor is not yet large enough to completely compensate for the smallness of m^Va as a diluter 
mass in Eq. (3.11), we show below that when both A''2 and are introduced into the game, 
the dilution factor S can be large enough for the correct DM relic abundance. Depending 
on the temperature of the QCD phase transition, there is a window of M^y^j for this effect 
to be significant. This relatively low scale window will be quantified in the next section to 
lie around 5 TeV or so; for the time being the essential point is just its existence. It is only 
a window, since for even smaller Mw^, below ~3 TeV, the freeze out temperature is below 
N2 mass, and the dilution factor gets Boltzmann suppressed. 

3.6 Diluters and dilutees: summary 

This is an appropriate moment to summarize our findings in a single spot and take a look 
at the parameter space of the LRSM in view of the relic density of the RH neutrino A^i, 
which we show in Fig. 3. The shaded belt regions show the possible parameter space where 
significant entropy production has a chance to take place. They correspond to the lifetime 
of the diluters within the 0.5-2 second range. If they are much heavier than GeV, they 
decay too fast, unless Wr is heavy (> 20 TeV). In this regime, the parameter space has a 
simple scaling law mj\i (x Notice that only when mjy ~ fn-w + lower values of Myy^ 
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Figure 3. Parameter space relevant for the warm dark matter in the minimal LR model. Left. The 
shaded regions (green, blue and red) labeled e, /i and r depict the regions, where the lifetime of the 
diluting N lies between 0.5 and 2 seconds, with N coupling predominantly to a single flavor. The 
spikes due to the phase space suppression in the decay of N allow for lower values of Mwr- Right. 
A zoom to the region Mw^ ^ lOTeV of our primary interest. Also shown are the theoretical lower 
limit on Mw,i from kaon mixing (which coincides with the lower limit set by the current LHC direct 
search of W e/n + missing energy), as well as the 14 TeV LHC reach. 

are allowed, producing the spikes seen in Fig. 3 in accord with the flavor diagonal in 
Eq. (3.15). 

Furthermore, the above regions terminate at values of for which the diluters start 
to feel the non-relativistic suppression during its freeze out. In the case of i = t, this gives a 
lower bound Afw^ > 16 TeV [31]. Once again, this is the reason why we focus here on lighter 
diluters which do not couple to r in order to have LR symmetry near the TeV scale. As 
we showed above, and as will be discussed in detail in the next technical session, the main 
message here is the existence of a window of the LR scale, potentially accessible to the LHC, 
besides the expected high scale scenario discussed in the past. 

Within the region below the belts, the decays of diluters are usually too fast and one 
ends up having too much DM in the universe. For the region above, however, the decays 
are not necessarily slow, because there is still a possibility to decay via the Dirac Yukawa 
couplings, which we have neglected so far. These couplings induce the mixing between the 
RH and SM neutrinos, and new decay channels — t- 3i^ or i'e~^e~ mediated by the SM gauge 
bosons open up. If these decay channels dominate over the SU{2)ji gauge interactions, the 
production and dilution of A^i crosses over smoothly to the z^MSM [30, 32] case in the phase 
diagram. 

However, there is a subtle difference from the usual z^MSM picture due to the presence 
of new gauge interactions, which now enter the game. Even though they do not play any role 
in the decay of the diluter, they may still thermalize and over-produce the DM candidate A^i, 
if the universe started out at a high enough temperature. In this case, a late-time dilution is 
still necessary, which calls for niN > 1 - 2 GeV [30] . This corresponds to the magenta shaded 
region in Fig. 3. Another way out would be to consider the reheating temperature of the 
universe after inflation to be sufficiently lower than the typical freeze-out temperature on the 
order 100 MeV - GeV. In this case, A''! has to be produced in a non-thermal way [27, 72, 73]. 
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4 The Boltzmann Approach 



As promised, we now come to the section of the technical aspects of our computations. 
Here, we implement the picture described in the previous section and numerically solve for 
the dilution of DM relic abundance. This quantifies the window of low scale LR symmetry 
consistent with a warm DM candidate. If not yet done, this is the right moment for our 
reader to take a (short) coffee break. 

In order to set the stage, let us recapitulate the main points which follow from the 
qualitative study of the previous section. 

1. The lightest neutrino A^i, the DM candidate, weighs around a keV, and therefore is 
always relativistic throughout the thermal history of interest here. It is coupled to 
r-lepton only {V^ = 1), which suppresses the charged current interactions when the 
temperature drops below the r mass, leaving room only for neutral currents mediated 
by Zlr. 

2. Such a flavor structure guarantees that a diluter does not decay to A^i if its mass is 
below rnr- 

3. There is a profound difference between the decoupling of the diluters N and the light 
neutrinos (and Ni). A diluter's freeze-out temperature is determined by the single 
scattering with quarks/pions through charge-current interactions. 

4. The late decay of a diluter happens at the temperature around MeV, which is well 
below their freeze-out temperature and therefore we are going to treat these processes 
in two separate stages. 

5. In order to optimize the entropy production, the diluters with lifetimes as long as 
a second should be as heavy as possible. In turn, this implies ~ ^ 1 and 
TTiATj ~ "T-TT + iTT'fijfnNs ~ mn + i^^e- In this case, the pionic decay dominates and the 

dependence is shown in Fig. 3. 

4.1 Freeze out 

We now turn to the study of the Boltzmann equations governing the freeze-out of RH neu- 
trinos. In order to keep track of the expansion of the universe, we define an arbitrary 
temperature T, which simply scales as 1 /R. One can think of T as the temperature of some 
fictitious relativistic species which freezes-out at some initial temperature Tj. Notice that the 
photon temperature is no longer a convenient choice as the number of relativistic degrees of 
freedom, i.e. f?* changes dramatically with time, especially during the QCD phase transition. 
In any case, one can always solve for the photon temperature Ty in terms of T, using entropy 
conservation 



In practice, we choose the initial temperature to be Tj = 10 GeV and find T < Ty < 2 T 
throughout the freeze out process of A^'s. 

The set of Boltzmann equations describing the freeze out of Aj are 



(4.1) 
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Figure 4. Solid lines representing the yields of Ni (magenta), N2 (blue) and (green) during the 
freeze out process, as solved using the Boltzmann equations in Eq. (4.2). In both panels, the masses 
of ^^2,3 correspond to the lifetime of 1.5 sec. Left. Yields of Ni as a function of temperature during 
the freeze-out process for a fixed value of = 5.5 TeV. The dashed curves corresponds to the 

thermal equilibrium yields of Ni, while the kinks reflect the reheating effect of QCD phase transition. 
Right. The reader may wish to take a deep breath before staring at the this panel. Final yields for 
different LR scales as a function of M\y^ computed using different values of Tqcd as explained in the 
legend (cf. with the naive estimate in Fig. 2). 



where z = uin^/T, = niN^/Ty and the entropy density and the Hubble parameter are 
defined as usual 

' = 1^5*^^7' H = 1-66^5:^ ■ (4.3) 

The single- and pair-annihilation reaction rates T/v^^jT^^^^^ and 7^^^. are given in the 
Appendix B. As discussed in the section above, the single-A^ processes mediated by Wr 
dominate by far over the pair annihilation processes. The latter interactions only affect the 
result by less than 5%. For LR symmetry scale close to TeV, the charged-current interaction 
rate of A^i is relatively smaller compared to that of A''2,3, due to the Boltzmann suppression 
in the r-lepton number density. 

We have numerically solved the above Boltzmann equations to calculate the yields of the 
three RH neutrinos. We start evolving the equations from a sufficiently high temperature, 
where all the states are in thermal equilibrium. A sample solution of the yields depending 
on z (temperature) during freeze-out is shown in the left panel of Fig. 4, where we take 
^Wr = 5.5 TeV, Tqcd = 350 MeV and we obtain masses of by fixing their lifetime 
equal to 1.5 sec. In this case, A^i freezes out just before the QCD phase transition and its 
temperature barely receives the heating from change. In contrast, N2 and N3 can freeze 
out after the transition and before they become non-relativistic, thus their number densities 
get enhanced. This realizes the large relative ratio discussed in Eq. (3.19). 

We proceed to compute the dependence of the final yields for a varying LR scale and 
show Y/v- as a function of in the right panel of Fig. 4. In order to obtain the final 

yields, we take the numerical solution at z ~ 30, where the freeze-out is basically finished, 
and where the time is still early enough so that have not yet started to decay. The final 

^We choose as tlie normalization point, since N2, being heavier than is the dominant source of 

dilution. 
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yields obtained here will serve as initial conditions for the Boltzmann equations describing 
the late decay of in the following section. 

The realistic numerical results of this section are to be compared with the estimate in 
Fig. 2. One can see that the largest ratios of Yiq^^/Yj^^ can be achieved for My/j^ between 
4 — 6 TeV. During the freeze-out process, even though the interaction rate has dropped below 
the Hubble rate, its effect on the evolution does not disappear immediately. This is why the 
final ratio tends to be smeared compared to the rough estimate in Fig. 2. 

4.2 Late decay 

The calculation of entropy production due to the late decay of diluters is based on the 
assumption that the microscopic interactions of the remaining SM particles in the thermal 
plasma are still fast enough. This guarantees that the relaxation time for the decay products 
to equilibrate with the plasma is much shorter than the time scale of the universe expansion. 
Therefore, the heat release dQ from the decay is completely transferred into the energy of 
radiation. The final entropy release A5 = J dQ/T [28] can be obtained by studying the 
Boltzmann equations for the energy density evolution [86]. In our case, the typical range of 
temperatures for the matter domination and decay to happen is between 0.5 — 10 MeV. The 
set of relevant Boltzmann equations for the energy density of matter and radiation are: 

dpR 



dt 
dpNi 

dt 
dpN2 

dt 
dpNi 



+ AH PR = T2P2 + Fgps , (4.4) 

+ 4HpN, = , (4.5) 

+ 3HpN, = -T2P2 , (4.6) 

+ 3HpN, = -FsPs ■ (4.7) 



dt 

The initial conditions are obtained by matching this late decay regime to the freeze-out 
regime in the previous section, at time 

(4.8) 
(4.9) 
(4.10) 
(4.11) 

We choose a matching point with = 30, or photon temperature around 15 MeV, 
where the temporary matter domination of diluters has not yet begun. The corresponding 
Hubble time can be calculated using t^ = 1/2H{T^^^). We evolve the above equations up to 
a time tgn when the entropy production has finished, typically much larger than the lifetime 
of a diluter. 

The reheating temperature after the decay is defined as the temperature when ra- 
diation again starts to dominate the universe, which happens around t = tm- It can be 
extracted from the radiation energy density at that time 

Pr{tn) = '^9*{Tr)T^ ■ (4.12) 
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Figure 5. Time evolution of energy densities for different species in the universe. We have used the 
same set of parameters as the left panel of Fig. 4. The universe temporarily enters a matter-dominated 
phase before A^2,3 decay. 



The dilution factor S is defined as the increase of total entropy S due to the late decay, 
as was discussed at length in Section 3. It can be derived based on the fact that radiation- 
like SM particles absorb all the heat from A^2,3 decay, therefore the products of the diluters 
dominate by far the entropy density. On the other hand, our DM A^i is by now completely 
decoupled from the plasma and since it is not a product of the heavier neutrinos, it simply 
dilutes as pjsf-^ oc R~'^. One then gets an improved dilution factor, which was estimated before 
in Section 3 in the sudden decay approximation (see Eq. (3.10)) 





_ s{tf) V{tf) _ 




3/4 




3/4 


'S'(im) 




_PR{tia) _ 









(4.13) 



This dilution factor will be used to rescale the relic abundance of A'^i in Eq. (3.6), where Y^r^ 
is also calculated numerically from the previous subsection. 

A sample solution with model parameters M^y^ = 5.5 TeV, Tqcd = 350 MeV and 
ttv = 1-5 sec is shown in Fig. 5. The initial conditions are obtained by matching to the sample 
values obtained from the freeze-out. We find the energy density of A'^2,3 can temporarily 
dominate over that of radiation before their decay, and the universe temporarily enters matter 
dominated expansion. The matter domination factor (pjy^ + PN-i)/ PR can be as large as 3. 
In this case, the reheating temperature is Tr ~ 0.7 MeV and could be consistent with helium 
abundance in standard BBN [87]. We comment on this in more detail, as well as on the 
implications for the CMB in the following section. 

In order to quantify the low-scale Mw^ window in which the relic density of A^i agrees 
with that of dark matter, we do a complete numerical study and compute relic density 
after the maximal dilution, with tn = 1.5 sec or ~ 0.7 MeV. In Fig. 6, f^TVi is plotted 
in together with the WMAP favored value f^oM (Eq. (3.7)) shown in green bands. For 
a varying lifetime, the reheating temperature scales as I/^/tn, while the dilution factor S 
scales as ^/TN■ We find that the desired window exists for mi^^ = 0.5 keV. In the case of nearly 
first-order QCD phase transition happening at 350 MeV, the light Wr window lies between 
4 — 8 TeV, while for a second-order transition at around 400 MeV, the window shrinks and 
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Figure 6. The minimum DM relic density that can be accommodated in the LR theory is plotted 
for two fixed DM masses as a function of Mwa- Solid, dashed and dot-dashed lines correspond to 
2^QCD = 350 MeV (nearly first order) and Tqcd — 150,400MeV (second order), respectively and a 
fixed lifetime t^^ ^ = 1.5 sec. The green bands from dark to light correspond to DM relic abundance 
at 1, 2, 3(7 confidence level from WMAP fit. Keep in mind that the DM relic density is linearly 
proportional to niNi- 

shifts to around 8 — 9 TeV. For lower QCD transition temperatures, the window moves to 
lower Miy^, but calls for an even lighter DM mass rrajv^ ~ 0.4 keV. For such a value, the 
window is there, regardless of the nature of the QCD phase transition. 

4.3 Dark Matter and LR scale: summary 

Here we summarize our quantitative results after numerically scanning over the parameter 
space. In Fig. 6 we plot the minimal relic density of the DM candidate A^i as a function 
of the LR symmetry scale M\y^ for two different values of m^^^ together with the allowed 
values coming from observations. Let us explain the qualitative aspects of the plot on the 
left panel of the figure, for the more conservative value rriN^ = 0.5 keV. 

First of all, the reader should recall that it was crucial to decouple our DM A^i as early as 
possible in order to separate its freeze-out temperature from the one of the diluters, profiting 
from the QCD phase transition. This is why A'^i had to be coupled to r in order to have 
predominantly weaker neutral interactions at that point in thermal history. In order for this 
to work, Wr can not be too light; otherwise the Tfs of all A^'s are equal and below Tqcd 
and A'^i remains over-abundant. 

Now for the window around My/^ w 5 TeV obtained by successfully separating Tf of the 
diluters and the dilutee. As seen in the figure, this works only for the nearly first-order QCD 
phase transition happening at around 350 MeV. If the DM mass is allowed to be this low, it 
is remarkable that the Wr mass window lies at the heart of the 14 TeV LHC reach [76, 77]. 

Further increase in the Wr mass results in the merging of these two different temper- 
atures, now above Tqcd- As we go along, taking bigger values of Wr mass, the relic DM 
abundance keeps falling down due to the fact that the mass of the diluters, directly controlling 
the dilution, gets increased, see Fig. 3. 

The DM relic density simply scales linearly with m^r^, therefore if one were allowed to 
go to even smaller values of rriTVi, say 0.4 keV as in the right panel, the correct DM abundance 
could be achieved for basically any value of Tqcd • In such case the Wr mass window either 
expands or the upper boundary completely disappears. 
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Constraints 








Dwarf Galaxy 


> 0.4-0.5keV 






Lyman-a 


> 0.5-lkeV 






BBN & CMB 




< 1.5 sec 




0u2p 






> 6-8TeV 


LHC-14 reach 


0-Mwn 




< 6.3 TeV 


A sample point 


0.5 keV 


1.5 sec 


4-7TeV 



Table 2. Various constraints on the masses and liftetimcs of relevant states within the LRSM, coming 
from astrophysical, cosmological and terrestrial experiments, together with a sample point in the DM 
scenario. 

For the sake of completeness, we wish to remind the reader that there was also the 
other possibility of A^i being coupled to light leptons, but this would force Wr to be heavy, 
definitely above the LHC reach and therefore not of our interest here. 

5 Further Constraints and the Uncertainties 

We now turn to astrophysical, cosmological and low energy constraints on the DM scenario 
discussed above. They can be classified into constraints on the DM mass mjv-^, lifetime of 
diluters tn and the left-right symmetry scale Mw^^. These bounds are summarized in Table 2. 

5.1 Dwarf spheroidal galaxies and Lyman-a forest 

Cosmological observations put a lower limit on the mass of the DM candidate. The most 
reliable and conservative bound comes from the study of dwarf spheroidal galaxies, whose 
content is commonly believed to be dominated by dark matter. If fermionic DM inside such 
astrophysical objects can be regarded as a degenerate Fermi gas, the requirement that the 
DM velocity on the Fermi surface must be less than the escape velocity, gives a lower bound 
on its mass rriDM > 0-468]'jgQg2 keV [65]. A more sophisticated analysis which compares 
the maximum phase space density [64-66] with observations gives a slightly stronger bound, 
"T-DM > 0-557^QQgY keV. Our case with mAr^ 0.5 keV is consistent with these lower limits. 

On the other hand, a lower limit on the DM mass can also be inferred from studying the 
absorption lines in the Lyman-a forest, which mainly constrains the maximal free-streaming 
length of warm DM. Most N-body simulations are carried out for structures formed at red- 
shifts approaching the non-linear growth regime, and are subjected to large uncertainties. 
The most recent constraints are made for warm RH neutrinos produced via non-resonant 
production, with Tjv^ ~ T^, with a range between 8-14keV [78, 79]. In our scenario, A^i 
freezes out earlier than SM neutrinos and gets further diluted. In such case, its temper- 
ature, free-streaming length and the mass lower bound should be reduced by a factor of 
{g*{Tf^^)/{g^{Tfi)S)f^^. Therefore, the lower bound for our case is mDM ^ 0(1) keV. We 
also notice that a recent analysis [80-82] using a thermal relic warm DM with mass 0.75 keV 
finds consistency with WMAP and Lyman-a observations. In view of the uncertainties men- 
tioned above, we stick to our conservative lower limit about 0.5 keV. 

5.2 CMB and BBN: neutrino thermalization 

Another class of cosmological constraints is related to the reheating temperature after the 
late decay of As we learned in this study, a large enough dilution implies a longer 
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lifetime of N and therefore a lower reheating temperature. If this temperature drops below 
the SM neutrino decoupling temperature, which is around 1 — 2 MeV, the thermalization of 
neutrinos becomes inefficient. In turn, this has a strong impact on both, the CMB power 
spectrum and the production of light elements during BBN. Naively, one could conclude 
that are almost no light thermalized neutrino species, but of course the diluters' decays into 
neutrinos lead to nonzero effective neutrino number as discussed below. 

CMB. The moment of matter-radiation equality determines the number of relativistic de- 
grees of freedom, which can be measured by observing the CMB power spectrum and is 
parametrized by the effective number of neutrino species, N^s. The determination of A^efT 
has been improved over the years with the most recent best fit at N^s = 4.34^Qgg (68% CL) 
from WMAP-7 [88] , while another recent analysis, combining the low redshift data with the 
5-year WMAP analysis, comes up with a slightly different value of N^g = 3.77 ± 0.067 (68% 
CL) [89]. 

In our scenario, both diluters decay into neutrino-rich final states. If it were not 
for these decays, the light neutrinos would have been extinct. Once again, this testifies a 
completely different history of the early universe from the standard one. More precisely, 
the lighter state decays as N3 — )• vr+e" — )• e"^e~f^z/^fe and the heavier A''2 — )• vr^/i" — )• 
e~ v^Ufj_i>^U(,i>e (or their anti-particles), where the average energy of the final state neutrinos 
ranges from 10 — 50 MeV. As a rough estimate we simply count the number of neutrinos 
produced in the decay. Using Eq. (3.5) with g*{Tf) ~ 20, appropriate for the freeze-out of 
A^2,3 below Tqc£), it implies that each of the A^2,3 number density counts as a half of the 
usual light neutrino one (recall that g^:(MeV) ~ 10). Since on average they decay into four 
light neutrinos, one gets roughly A'efr — 4, however, we expect that a more precise numerical 
solution would give a somewhat lower value. Namely, the weak interaction cross-section of 
neutrinos scattering on electrons or protons is enhanced by the neutrino energy [90] , compared 
to the thermal one, and therefore these neutrinos down-scatter until they lose most of their 
energy to the plasma [91]. They will also start to annihilate with each other when becoming 
sufficiently populated [92]. A complete analysis of neutrino thermalization in this scenario is 
beyond the scope of this paper. 

BBN. A late decaying particle resulting in a low reheating temperature before the onset 
of BBN, could drastically change the prediction for the primordial Helium abundance by 
affecting the neutron-proton number ratio. 

The injection of energetic pions resulting from the decay would increase the neutron- 
proton ratio. Notice however, that the pions in this scenario are rather soft, since misf lies 
near the threshold, therefore they do not scatter with nucleons before decaying [94] and the 
above problem can easily be evaded. Second, the lack of neutrino thermalization changes 
both, the contribution of neutrinos to the Hubble rate, and the average weak interaction rate 
with competing effects. The analysis of [87, 92, 93] shows that a reheating temperature as 
low as 0.7 MeV could still be compatible with the observed Helium abundance. 

5.3 Neutrino-less double-/? decay 

As mentioned in the introduction, the Majorana nature of the RH neutrinos and the as- 
sociated lepton number violation (LNV) plays a crucial role in the study of dark matter 
in the minimal LR model. The textbook example of LNV is the neutrinoless double beta 
decay (Oz^2/3), often erroneously associated with neutrino mass only in spite of the fact that 
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more than 50 years ago it was argued that new physics may be equally responsible for this 
process [95, 96]. 

The LR theory is tailor-made for the new physics contribution to the Oi^2/3. Since left 
implies right, there must of course be a right-handed counterpart to the usual left-handed 
neutrino contribution [35]. Recently, an in-depth study has been performed [14], which 
emphasized the profound connection between Oz^2/3, LNV at colliders [10] and lepton flavor 
violation (see also [97] for a recent study of Oi^2/3 in LRSM and [98] for a study on future 
linear collider signals). This leads to a derivation of a serious constraint [15] on the scales of 
LR symmetry, which depends on the flavor structure of V^. In our case, since ends up 
being diagonal, the limit becomes quite acute. From Fig. 1 in [15] one can conclude that for 
the preferred value of ttt-tv around 140 MeV, the mass of Wr should lie above ~6-8TeV. 
Given the large uncertainties of the nuclear matrix elements [99, 100], we cannot claim with 
sufficient certainty that this setup is completely outside the range of the LHC and we leave 
the precise determination of this particular bound to the experts. In any case, upcoming 
experiments are about to probe this region [101, 102]. 

6 Conclusions and Outlook 

The final outcome of this study is quite simple and striking. The minimal left-right symmetric 
model can naturally provide a warm dark matter candidate in the form of a keV right-handed 
neutrino, with a hope of observing the RH charged gauge boson at the LHC. For those 
of you who do not care about the nitty-gritty of the theory and our analysis, this is the main 
message we would like you to take home. 

For this to be true, the following spectrum emerges, rriNi — keV (Dark Matter), m]\j^ 
rriTT -|- m^, rriN.^ + (Diluters); together with a particular flavor-diagonal structure 

yli yR^ yR i This is the only option to have a reasonably light Wr accessible at 
the LHC. We find on top, for mjy-^ lying around 0.5 keV, a narrow window of around 
5 TeV, otherwise has to be larger than about 20 TeV. We find this isolated window very 
interesting. 

A noteworthy fact. For this DM picture to work, with the spectrum taken as above, 
at no point we need to assume that the universe ever reached very high temperature. The 
highest temperature we talk about here is less than about GeV, a modest extrapolation of 
the BBN temperatures. This puts the warm DM picture of the Left-Right Theory on quite 
firm grounds. 

The million dollar question is: how to test this DM scenario of the LRSM? Ideally, 
one would like to directly measure the masses and mixings of the right-handed neutrinos 
together with a mass of W/j. The only way to measure uin and at the LHC is through 
the golden KS channel [10] of two charged leptons and two jets without any missing energy. 
Unfortunately, this requires much larger masses mjv ^ 10 GeV than the resulting spectrum 
above, which could manifest itself simply as the missing energy due to the extremely long-lived 
diluters. At this point, it seems difficult to imagine a conceivable way of directly measuring 

and/or niN, despite the fact that the resulting parameter space in the right-handed 
leptonic sector is well determined. On the other hand, one can search for indirect signals 
using low-energy processes, similarly to the sterile neutrino case [103-105] (see also [106- 
109]), which have a correlated dependence on both, the mixing parameters and the mass 
spectrum. In view of the smallness of the Dirac Yukawa couplings, the presence of new 
heavy gauge bosons is more than welcome for the sake of visibility of these elusive, otherwise 
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sterile neutrinos. Moreover, a Wpt in the ~ 5 TeV window can well be probed at 14 TeV 
LHC in the case when light RH neutrinos escape detection with a potential to establish the 
chirality of the outgoing fermions [76]. ^ 

On a positive note, we predict a sizable rate for the neutrino-less double beta decay, 
already on the edge of exclusion by the current data. It is perhaps even more difficult to 
imagine that Oz^2/3 decay would not be seen in the currently on-going and planned experi- 
ments, if Wr were seen at the LHC. In this sense, our work provides an additional impetus 
for the dedicated Oz^2/3 decay searches. Moreover, although it is hard to verify precisely the 
picture that emerges from our analysis, it is easy to kill it. The smallness of RH neutrino 
masses and the absence of leptonic mixing, eliminates completely lepton number violation 
at colliders and any lepton number violation. Observing any such process would invalidate 
completely the scenario offered here and dark matter would need another origin. 

In our opinion, there are still some issues that deserve deeper understanding and more 
quantitative study. The foremost is the neutrino thermalization before the BBN epoch, which 
affects both the light element production and the effective neutrino species A'eff measured 
by CMB (see e.g. [91]). This could seal the fate of our scenario, especially in view of better 
sensitivity of the upper coming data [110-112]. Furthermore, there seems to be a lower limit 
on the X-ray line flux from the decay of our warm dark matter, due to the flavor structure 
imposed by left-right symmetry and the constraints from our dark matter analysis. 
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A Annihilation Cross Sections 



In order to write down the Boltzmann equations for the freeze-out process, we summarize the 
main processes governing the freeze-out. In order to get the thermal reaction rates discussed 
in Appendix B, one needs to know annihilation cross sections. We work with the flavor 
structure where A^i mainly couples to r. When the temperature drops well below m,-, all 
charge current interactions of A'^i become negligible due to Boltzmann suppression. 

One actually uses the so-called reduced cross section (we consider 2 — 2 processes 
1,2 — )• a,b only), defined as 



s J (2^2K (^^^'"^ ^ Hpi+P2-Pa-Pb)\M\-.{A.l) 



a) Single-N annihilations, Wji-exchange 

The single- annihilation processes are mediated by Wji in both s- and t-channels. In 



*Still, one of the authors of the paper (GS) is deeply disturbed by the sad outcome of this work, which 
says that the golden KS channel will not be measured at the LHC. 
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the case massive A''2,3, the reduced cross sections are 

a, {Ndn ^ enun){s) = ^^^^^ , (A.2) 

t \ H H H)\ } 487rM4,^ ' ^ ' 

c^s "-{NeR eRdR){s) = . , (A.4) 

where qr is the SU{2)r gauge coupling, A'^c = 3 is the color factor and s is the center- 
of-mass energy squared of the 2—7-2 scattering process. Our results agree with those 
calculated in Ref. [113], up to the number of generations. Since we are discussing physics 
with energy scale around and below GeV, the 1/M^^ expansion has been made in the 
above expressions. In the Boltzmann equations, we also take into account of the charge 
conjugation processes. 

In the case of Ni, which is effectively massless, we have to keep the mass of r lepton 
to which it predominantly couples. The corresponding cross sections can be obtained by 
simply making the replacement mi\f — )• mr in Eqs. (A.2)-(A.4). 

b) Pair-of-N annihilation to the same flavor charge leptons, Wr and Z^R-exchanges 

The pair-annihilation process can take place through t- and u-channel Wr exchange, as 
well as s-channel Zlr exchange. For the case of A''2,3, the final leptons are massless and 
the cross section is 



g {3 cos 26w — 1) sec 9w sec 20w 



1927rM2 m2 



_^ g'^ sec^ 9w 



- 4m^)^/^ (^s - 4m|r + (s + 2m%) cos 26*1^^ 
16 cos^ 20w{2 cos 2ew - 1) ((« - m% f - + 



30727r(l - tan2 ew)^M^^^s 

+{8cos26w — 5cos40vF — 7)(m^ — s) ^cos4^vf(s + 2m^)(2s + m%) + 6m% — 9m%s + Qs^ 

(A.5) 

For the case of A'^i, where the mass of r in the final states cannot be neglected. Here the 
t-, u-channel interference is suppressed by the smallness of A'^i mass. The cross section is 



967rM2 Ml 

Wr ■^lr 



^/s(s - 4m^)^/2 {3s - 6ml - {s - 4ml) sec26'iy) 



+ 



g^ sec^ Owifnl — s) 



12(12 cos 2ew + 5(2 cos 2ew + cos AOw + 2) sec^ Ow 



7687r(l -tan2 ew)^M^^^s 

36)mls + (8 cos 26'vf - 5 cos A0w - 7) sec^ (9vf(1 + tan*^ Ow) {2s^ - mis + 2mf) 

(A.6) 
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c) Pair-of-N annihilations to other fermions, Z^fi-exchange only. The cross section of the 
pair-annihilation process can also take place through s-channel Zlr exchange 



a^'^.^iNN rj%s) ^^^^ 



48 cos^ 261^/(8 COS - 1) 



91367r(l - tan2 OwYM^^^s 
X {{s - m%f - + - (32cos26lvy + 53cos46ivi/ + 63)(s-m^) (A.7) 

X ^cos46'vk('S + 2m^)(2s + m\) + — 9m^s + 6m^ 

where /' is any fermion except for a charged lepton that couples to and Wn (recall that 
the RH leptonic mixing is flavor diagonal, as discussed throughout the section sec. 3). For 
the case of A'^i, one can simply take the limit m^v — >• 0. 

B Effective Thermal Rates 

The processes involving the inter-flavor pair annihilation of A^'s via Wr are always subdom- 
inant and have been neglected. For 2—7-2 processes, the reaction rate weighted by thermal 
distribution is 

7(1, 2^ a, 6)^-^ r dsa{s)^sKJ^\ . (B.l) 



647r4 . ^ / 

Here, the photon temperature is used because the thermal averaged rates are for those 
particles within equilibrium with the photon and s = Min {(m-i + 7712)^, (m-a + ^b)"^}- 

We name the reactions rates J^'^ for the sum of processes in a), for process 

b) and 7^^jv process c) of the Appendix A. Here, we list the thermal reaction rates for 
single- A^ annihilations. 

^""^ ~ 16^5m4,^4 ' (B.2) 

where = 171^2 /T^ and the pre-factors 2 represent the contribution of charge conjugation 
processes. 

In Fig. 7, we plot the the thermal averaged reaction rates for the processes that control 
the freeze out of A^i and N2, for Mp^/^ = 5TeV, m]\f^ = 0.25 GeV. For the case of A'2,3, 
the single- A^ processes (with rate 1^^^) always dominate. On the other hand, for A^i, the 
single- A^ interaction receives the Boltzmann suppression due to the presence of r lepton and 
become subdominant for temperature below 200 MeV. 

In Fig. 1, we plot the the thermal averaged reaction rates that dominate the freeze out 
of A^i and A^2 , for the same set of parameters as in Fig. 7. We also plot the Hubble expansion 
rate multiplied by the thermal number density of A^2- Using the naive decoupling condition, 
7 = nH, we find the freeze out temperatures are Tyi ~ 450 MeV and Tf2 ~ 250 MeV, 
respectively. This difference makes the scenario, where A^i freezes out before and A^2 after 
the QCD phase transition (the tilde on the nH curve), possible. 
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1 GeV/r^ 1 GeV/r^ 




1 GeV/r^ 

Figure 7. Thermal averaged reaction rates for the processes that controls the decoupling of Ni (red), 
7V2 (blue) and (green), for parameters Mwrc — 5TeV, rriM^ = 0.25 GeV, niN^ = 0.14 GeV. For 
each N, the relative value of different rates does not depend on the choice of M^y^. 
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